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Abstract Electromagnetic suspension system (EMS) is mostly used in the field of high-speed vehicle. In this paper, 
a space exploring vehicle quarter electromagnetic suspension system is modelled, designed and simulated using 
linear quadratic optimal control problem. Linear quadratic Gaussian and linear quadratic integral controllers are 
designed to improve the body travel of the vehicle using bump road profile. Comparison between the proposed 
controllers is done and a promising simulation result have been analyzed. 
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1. Introduction 

Electromagnetic suspension (EMS) is the 
magnetic levitation of an object achieved by using 
continuously altering the power of a magnetic field 
produced through electromagnets the usage of a 
remarks loop. In maximum cases the levitation impact 
is commonly due to permanent magnets as they don't 
have any strength dissipation, with electromagnets 
simplest used to stabilize the impact. According to 
Earns haw’s Theorem a paramagnetic ally magnetized 
body can't relaxation in solid equilibrium whilst 
located in any aggregate of gravitational and magneto 
static fields. In those forms of fields an equilibrium 
circumstance exists. Although static fields cannot 
provide stability, EMS works by way of always 
altering the modem-day sent to electromagnets to 
change the energy of the magnetic area and permits a 


solid levitation to occur. In EMS a feedback loop 
which constantly adjusts one or extra electromagnets 
to accurate the item's movement is used to cancel the 
instability. Many structures use magnetic attraction 
pulling upwards towards gravity for these styles of 
structures as this gives some inherent lateral stability, 
however a few use an aggregate of magnetic appeal 
and magnetic repulsion to push upwards. Magnetic 
levitation generation is important as it reduces 
electricity consumption, in large part obviating 
friction. It additionally avoids wear and has very low 
protection necessities. 

2. Mathematical Modelling of the 
Electromagnetic Suspension System 

Figure 1 shows a quarter vehicle electromagnetic 
suspension system and the electromagnetic tire design. 
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Figure 1 Quarter vehicle electromagnetic suspension system and the electromagnetic tire design. 
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The tire is suspended with an initial length xO for an 
initial current iO. The electromagnet tire has a 
potentiometer fixed inside it and at the layer of the tire 
there is a coil connected. The potentiometer has a 
source voltage E and a resistor length D. When a road 
disturbance is accrued, the tire with the attached 
potentiometer rod will move upward and downward 
with a length y so the metallic body mass will move. 
The design is as follows. 

Apply Kirchhoffs voltage equation for the 
electric circuit 


v = v B + v r 


, , di 

u(t) — iR + L — 
v ’ dt 


( 1 ) 


Where u, I, R and L is applied voltage input, 
current in the electromagnet coil, coil’s resistance and 
coil’s inductance respectively. 

Energy stored in the inductor can be written as 


^I Stored ~ 2^ 


( 2 ) 


Since power in electrical system (Pe) = Power in 
the mechanical system (Pm), 

dW, a 


Where 


P =- 


dt 


f 

J m 


dW, 


I Store 


dt 


p„, 

dx 

and 

dt 

dt 

dW IStore 

dx 

dx 


( 3 ) 


Where fm is known as electromagnet force Now 
substituting (2) in the equation (3), 


-Li 2 

2 


_ d 

dm 7 

ax 


2 dx X 


( 4 ) 


Since the inductance L is a nonlinear function of 
body travel position (x) we shall neglect the leakage 
flux and eddy cuii'ent effects (for simplicity), so that 
the inductance varies with the inverse of body travel 
position as follows: 


L = — Where in K = 
x 2 


(5) 


Where, p.0 is the inductance constant, A is the 
pole area, N is the number of coil turns and k is 
electromagnet force constant. 


1 , 2 d f k 


^ m 2 dx 1 x 




2 


f ,-2 


( 6 ) 


If fm is electromagnetic force produced by input 
current, /g is the force due to gravity and / is net 
force acting on the vehicle body, the equation of force 
can be written as 


fg=fm +f 
= fm + m 


d~x 

~dd 


dv 


'■ mg- 


‘(j) 

x(t) 


( 7 ) 


Where m = vehicle mass and v =dx/dt =dh/dt, 
which is velocity of the vehicle body movement. 

At equilibrium the force due to gravity and the 
magnetic force are equal and oppose each other so that 
the vehicle body levitates, i.e. fg = -fm and /=0. On 
the basis of electro-mechanical modeling, the 
nonlinear model of magnetic levitation system can be 
described as follows: 

The general form of an affine system 

^=f{z)+g(z)M (8) 

Is obtained by denoting variables for state space 
representation as follows 


dx 


Z 2 = 


dt 


■ = v 


z 3 = z 


( 9 ) 


Substitute equation (9) or the state variables in to 
equation (1) and (7) 

u{t) -Z y R + L ,xz 3 


mz 2 = m.g - - 


K 


r \ 


\ z \ j 


( 10 ) 


Then the nonlinear state space model is 
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z, = x 1 


2 2 = \ 8 


2 m J 


f \ 2 
Z 3 


u R 
z, =- z, — 


( 11 ) 


L L 

The fm is electromagnetic force produced by 
input current is related to the current and the road 
disturbance displacement will be 

fm=g( X > i ) ( 12 ) 

Using Taylor series linearization technique we 

have 


/ = 

J m 


f Sg} 

K dx j 


x - 


v Si J 


2 ki n 


2 ki n 


i = -- 


Where S = g( x A aml (. x 0 ,/ 0 ) 

operating point’s initial inputs. 

The applied force to the mass become 




is the 


M 


d 1 x 

dt 2 


•fm 


2 kL 


-x + - 


2k i 


°i 


(14) 


Then the linearized state space model becomes 


(z ^ 


V Z 3 J 
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mXn 
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v z 3y 


( 13 ) 


From the potentiometer 

v = £— (15) 

D v ’ 

So the final state space model becomes 


z i 


\ Z U 
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The system parameters are shown in Table 1 
below 


Table 1 System parameters 


No 

Parameter 

Symbol 

Value 

1 

Mass of the vehicle 

m 

1 kg 

2 

Coil resistance 

R 

G 

O 

3 

Coil inductance 

L 

0.2 H 

4 

Initial current 

h 

0.8 A 

5 

Initial displacement 

x 0 

0.03 m 

6 

Electromagnet Constant 

k 

Nm 

A 1 

2.9xlO A -6 A 

7 

Potentiometer voltage 

E 

5 V 

8 

Potentiometer distance 

D 

0.13 m 


The transfer function become 

VO _i_ 

Y(s) j J +50V+0.1375i + 6.874 

The state space representation becomes 
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3. Proposed Controllers Design 

3.1 LQG optimal controller Design 

LQG computes an optimal controller to stabilize 
the plant G (s) 


x = Ax+Bu +E, (16) 
y = Cx + Du + 6 

And minimize the quadratic cost function 


Jlqg lim j j[ x 11 ] 


T —ko 


'Q 

N T C R 


dt (17) 


The block diagram of a quarter vehicle 
electromagnetic suspension system with LQG 
controller is shown in Figure 2 below 


Plant Noise % (t) Measurement Noise 6 (f) 



Figure 2: Block diagram of a quarter vehicle electromagnetic suspension system with LQG controller 


The solution of the LQG problem is a 
combination of the solutions of Kalman filtering and 
hill state feedback problems based on the so-called 
separation principle. 


E 


[IT(<)‘ 

ihd. 


[£(0^O)] r j 


The input variables W and V are 


The plant noise and measurement noise are white 
and Gaussian with joint correlation function 





S(t-r) 


(18) 


W = 


Q N c 

n t c r 


;v = 


E N t 


N f 0 


The final negative-feedback controller becomes: 

F (&-= 


A—K f C 2 — B 2 K c + K f D 22 K c K f 


K r , 


0 


For the Gaussian noises ■ and 0 : 

H = 0.0005 and 0 = 0.0000001 

The value of Q and R is chosen as 

r 5 0 0" 

0 5 0 and R = 10 

V° 0 5 v 

The LQG controller becomes 


0.5678s 3 + 2.8754s 2 +1.43225 + 2,0123 
s 4 + 12s 3 + 35s 2 + 94s+ 112 

3.2 Linear Quadratic Integral Controller 
Design 

LQI computes an optimal state-feedback control 
law for the tracking loop. Block diagram of a quarter 
vehicle electromagnetic suspension system with LQI 
controller is shown in Figure 3 below. 


GcM 
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Figure 3 Block diagram of a quarter vehicle electromagnetic suspension system with LQI controller 
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Figure 4 Body travel response for a bump road disturbance 

For a plant sys with the state-space equations 

x = Ax + Bu (19) 

y = Cx + Du 

Tire state-feedback control is of the form 

u=-K[x,x^ (20) 

Where x t is the integrator output. This control 
law ensures that the output y tracks the reference 
command r. For MIMO systems, the number of 
integrators equals the dimension of the output y. LQI 
calculates the optimal gain matrix K, given a state- 
space model SYS for the plant and weighting 
matrices Q, R, N. 

The value of Q, R and N is chosen as 
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The LQI optimal gain matrix becomes 

AT = [0.2004 9.8905 1.0844 - 0.7071] 


4. Result and Discussion 

4.1 Body Travel Output Specification 

One of the major specification of a suspension 
system is whatever the road disturbance input, the best 


Body Travel Response for a Bump Road Disturbance 
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design performance of the body travel vertical 
displacement is to approach to zero. 

4.2 Comparison of a Quarter Vehicle 
Electromagnetic Suspension System with LQG and 
LQI Controllers for a Bump Road Disturbance. 

The quarter vehicle electromagnetic suspension 
system with the proposed controller’s comparison for 
a 10 cm bump road disturbance input simulation result 
for body travel response is shown in Figure 4. 

The simulation result numerical value is shown 
in Table 1 below. 


Table 2: Numerical values of the body travel 
simulation output 


No 

Systems 

Bump 

1 

Road Profile 

0.1 m 

2 

LQG 

0.01 m 

3 

LQI 

0.05 m 


Table 2 shows that the quarter vehicle 
electromagnetic suspension system with LQG 
controller body travel is minimum and improved the 
road handling criteria. 

5. Conclusion 

In this paper, a quarter vehicle electromagnetic 
suspension system design and analysis have been done 
using Matlab/Simulink. In order to increase the 
performance of the electromagnetic suspension 
system, LQG and LQI control technique is used. The 
main aim of this paper is to control the vehicle body 
travel based on road disturbance input. Comparison of 
the quarter vehicle electromagnetic suspension system 
with LQG and LQI controllers to improve the 
performance of the body travel output using a bump 


road profile. The quarter vehicle electromagnetic 

suspension system with LQG controller body travel is 

minimum and improved the road handling criteria. 
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